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Abstract 

We study the deconfining transition of SU(iV) gauge theories in 2+1 dimensions for 2 < 
N < 8. We confirm that the transition is second order for N < 3 and first order for N > 5. 
For the more delicate case of SU(4) all our evidence points to the transition being weakly first 
order. After extrapolating to the continuum limit, we obtain a deconfining temperature that 
can be well fitted by T c /y/a = 0.9026(23) + 0.880(43)/iV 2 for all N >2. 



1 Introduction 



In 2+1 dimensions, SU(iV) gauge theories have a coupling g 2 that has dimensions of mass. 
Thus the dimensionless coupling appropriate to physics on a length scale I is g 2 l: that is to say, 
the theory becomes free in the ultraviolet, and strongly coupled in the infrared. Moreover the 
theory is found to be linearly confining (see [1] and references therein) although no analytic 
proof exists. These properties are reminiscent of D = 3 + 1 gauge theories, and so D = 2 + 1 
gauge theories have received much attention in the hope that they might provide a useful 
stepping stone between the soluble D = 1 + 1 case, and the physically interesting D = 3 + 1 
case. 

In this paper we calculate some properties of the deconfining transition in D = 2 + 1 
SU(iV) gauge theories. We use standard lattice methods, extrapolate to the continuum limit, 
and then perform a large N extrapolation that enables us to extend our results to all values 
of N. These calculations parallel similar calculations in D = 3 + 1 [2]. 

There exist older calculations for SU(2) and SU(3) (see for example [3]) as well as more 
recent work on SU(4) [HE1E] and SU(5) [HE] most of which appeared while this work was in 
progress. Most of our calculations have appeared in [8] which contains various details omitted 
in this paper, and a summary of our early results was presented in [5]. The main change 
in this paper is that we use a more accurate calculation of the string tension pQ to provide 
the scale in which to express the quantities we calculate, and we examine the delicate case of 
SU(4) in greater detail. 



2 Lattice setup 

Our calculations are performed using standard lattice tachniques. We work on periodic cubic 
L 2 x L t lattices with lattice spacing a. The degrees of freedom are SU(iV) matrices, Ui, assigned 
to the links I of the lattice. The partition function is 

Z(ft) = j HdUie-f 3 ^-^^} (1) 

where U p is the ordered product of matrices around the boundary of the elementary square 
(plaquette) labelled by p. In the continuum limit 

2N 

P = — (2) 
ag 2 

where g 2 is the coupling, which in D = 2 + 1 has dimensions of mass. This is the standard 
Wilson plaquette action and the continuum limit is approached by tuning a — > and hence 
(3 = 2N/ag 2 — > oo. One expects that for large N physical masses will be proportional to the 
't Hooft coupling A = g 2 N [9], and this is indeed what one observes [ID] . So if we vary ft oc iV 2 
then we are keeping the lattice spacing a fixed in physical units, to leading order in N. 
When L s — > oo at fixed L t and fixed ft, the field theory is at a temperature 
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As is conventional, we use eqn([3]) to define T at finite L s ^> L t as well. By increasing (3 at 
fixed L t we increase T so as to locate the critical value f3 — {3 C at which the system deconfines: 
T c (a) = l/a(p c )L t . (We shall often write this as T c if there is no ambiguity.) An 'order' 
parameter we shall use is the plaquette, u p , averaged over the space time volume, 

up = jj- u p ; u v = ^ ReTrf/ P ( 4 ) 

p p 

where N p = 3L 2 L t is the total number of plaquettes. Another parameter is the Polyakov loop 

u 

l p (x) = Tr]JU t (x,t) ; l p = — ^l p (x). (5) 

t=l S X 



3 Methodology 

In this Section we briefly describe how we calculate the properties of the deconfining phase 
transition. For a more detailed review, with references, see for example [11] . 
Consider the specific heat on a L 2 s L t lattice, defined conventionally by 

j 2 C{{3) = §p(u p ) = N p (ul) - N p (u p f. (6) 

Suppose that there is a first order deconfining transition at (3 C so that (u p ) jumps from (u P;C ) 
t° (u p ,d) as we increase (3 through (3 C . On our finite spatial volume, V, the jump becomes 
a crossover. Let us choose the maximum of C{(3) as defining the critical (3 at finite volume: 
(3 c iy). Up to the weakly varying factor of (3~ 2 this coincides with the point at which (u p ) 
changes fastest with (3. It is easy to see that for large enough V this is the value of (3 at which 
the system is equally likely to be in the confined and deconfined phases and that 

lim ±- 2 C(f3 c (V)) = ^ «u PiC ) - (u p4 )f oc VLl (7) 

Thus at a first order transition the specific heat diverges linearly in V and the limiting value 
of C(/3)/V is directly related to the latent heat L h . Now, the crossover exists over a region of 
T where there is a finite probability for the system to be in both the confined and deconfined 
phases. Assuming that the free energies of the two phases are smooth as one passes through 
T c , this implies that the crossover extends over a range of 0(1/V) in T and hence that any 
definition of T c within this range will generically have a leading 0(1/V) correction at large 
V. So we expect that as V — » oo 

T c (oo) - T C (V) &(oo) - (3 C (V) h h (U 



T c (oo) &(oo) VT? ' \L 

(where we have used T = l/aL t = (3g 2 /2NL t ). Moreover, because the free energy oc iV 2 , we 
also expect that 

i N-*oo 1 

h K ( 9 ) 
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Thus if we calculate (3 C (V) from the maximum of C(/3), and do so for several V, we can use 
eqn(jHJ) to extrapolate to V = oo, thus obtaining f3 c (V = oo). We can then use the results 
of [1] to obtain the string tension at this /3, and hence T c (a) / y/cr(a) at a = l/L t T c . We now 
repeat the exercise for several other values of L t and perform an extrapolation of this ratio 
to the continuum limit. In the same way we can use eqnfJZJ) to obtain the continuum value of 
the latent heat, e.g. L h /Tj?. 

If the transition is second order, there is no discontinuity in (u p ) but there is a diverging 
correlation length. Rewriting eqn([6]) as an integrated correlation function 

j 2 C{P) = £<fe - (u P )) (u P0 - (u p ))) (10) 

V 

(where po is an arbitrary plaquette) we see that at the transition C((3) will diverge as V —* oo, 
with a power that depends on how the correlation length diverges with increasing V, i.e. on 
the critical exponents. The maximum of C (/3) is where the correlation length is largest (up to 
weak corrections) and provides a natural definition of (3 C (V). Using a conventional notation, 
the leading large volume behaviour is 

l^C(f3) = c V^+ Cl (11) 



and 

T c (oo)-T e (V) _ e (oo)-p c (V) _ h _ h (L t y d (u) 
T c (oo) &(oo) (VT^ \L S J 1 ] 

We can use this to obtain /3 c (oo) from our finite V calculations, then use [TJ to transform this 
into a value of T c (a) /y/a(a), then repeat all this for several values of L t , and finally extrapolate 
to the continuum limit. 

In practice using the plaquette and the specific heat is not the most accurate method. For 
our second order transitions the excitation whose mass vanishes typically has a weak overlap 
with the plaquette, so that it contributes a small peak to C{(3) unless the volume because 
extremely large. For our first order transitions the discontinuity is on physical length scales, 
so ((u Pt c) — (u P) d)) ~ 0(a 3 ) and once again one has to go to very large V to see a strong effect 
in C( S). 

For these reasons it is usual to also use the susceptibility of the Polyakov loop 

£ = <IM 2 > " <IM> 2 - (13) 

In the confined phase the Polyakov loop fluctuates about zero while in the deconfined phase it 
fluctuates about a non-zero value proportional to an element of the center of SU(iV). One uses 
\lp\ rather than l p since even in the deconfined phase the latter will (strictly speaking) average 
to zero for finite V through tunnelling between the N (Euclidean) vacua. While taking the 
modulus might seem like an ad hoc fix, in fact it is closely related to the adjoint Polyakov 
loop: TrA^p oc |Trp/ p | 2 — 1. The latter is a more natural and general 'order' parameter than 
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the usual fundamental Polyakov loop. For example, it will serve equally well if we decide to 
express the gauge fields in the adjoint reperesentation. For large N it vanishes as 1/N 2 in 
the confined phase, and is unsuppressed in the deconfined phase. For these reasons it would 
make sense to use a susceptibility based on the adjoint Polyakov loop. However, given that 
the definition in eqnf ll3l) is now customary, we shall adhere to it in this paper. Just as for 
C((3), we can define a (3 C (V) from the maximum of x/V and extrapolate it to V = oo using 
eqn([5I) or eqn (|T2|) depending on the order of the transition. The value of \jV at its maximum 
is not directly related to a physical quantity, unlike C((3), so we do not consider its infinite 
volume extrapolation. 

The reason that %jV provides a better signal for the phase transition than C{f3) is that 
the ground state of the flux loop that winds around the temporal torus becomes massless as 
T — ► T~ if the transition is second order. (This mass is an eigenstate of a transfer matrix 
along a spatial direction.) The Polyakov loop is an operator with a projection on this state. 
The fact that near T c excited flux loops are relatively massive, and that the projection on 
the ground state is quite good, for the range of a we work with, means that the effect of the 
vanishing mass on xjV is very pronounced. For a first order transition, on the other hand, the 
discontinuity in the average Polyakov loop is not directly related to a physical energy density 
and so is not as small as for the plaquette. As (3 — > oo Polyakov loops do get renormalised 
to zero, and using xjV becomes more difficult, but for our range a G [1/2T C , 1/5T C ] these 
problems are not severe and accurate calculations are easy to perform. 

We therefore base our calculations of (3 C on x/V- However we do use C{(3) to obtain the 
latent heat in those cases where the extrapolation in eqn([7|) can be usefully performed. To do 
so we recall that the energy density is given by 

I = -l^ aZ = Hp nZ , (14) 
V V d± V dT v ; 

where we use Z = ^ s exp{— i? s /T}, while 

d dT d dV d 



where the first equality in eqn fU5l) uses the expression in eqn([T]) for Z. Noting that the volume 
derivative gives the pressure, and that the pressure of the two phases is equal at T = T c , we 
see that the contribution of this second term to 

A Up = (u PjC ) - (u p4 ) = (u PjC ) - (up >d ) (16) 

will vanish at T = T c . Thus the latent heat, L^, which is the difference of the energy densities 
in the two phases at T = T c , can be written, using eqns fll4|15lfT6l) and T = l/aL t , as 

^ = -3L 3 t a^A Up ~ 3L 3 t (3A Up (17) 

where at the last step we use the asymptotic expression (3 = 2N/ag 2 . As we see from eqnfJTj), 
the quantity A u , and hence the latent heat, is given by the V — > oo limit of the specific heat. 
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The expression in eqn( U7l) is not unique. We can do better by noting that ad/3 /da = 
(afi)dj3/d(afi) where we choose for fi some physical quantity whose lattice value has been 
accurately calculated over a wide range of (3 (e.g. the string tension). Now parameterise the 
variation of a/i with (3 as a series in 1//?, evaluate (a/i)<9/?/<9(a/i) and substitute in eqn(TTTl). 
This will improve over using (3 = 2N/ag 2 by at least an additional 0(1/(3) correction. Having 
said that, we shall use the cruder estimate in eqn (fT7]) because our relatively large errors on 
L h do not demand a more sophisticated treatment. 

For some L t our range of V is not large enough to calculate from the specific heat. In 
those case we take a sequence of lattice fields generated close to T = T c , and identify long 
subsequences that are in either the confined or deconfined phase, using the value of l p to do 
so. We then obtain (u PtC ) and (u Pt d) from the appropriate subsequences, and hence A Up . This 
direct method of calculating A Up has the advantage of weak finite volume corrections and so 
does not require an explicit extrapolation to V = oo. However there are systematic errors to 
do with the assignment of lattice fields to a particular phase, and these errors only become 
small on large volumes where tunnellings between phases become relatively rare. Finally, 
using eqn ()17l) we transform our value of A Up into a value for the latent heat. 

An aside. Although one might imagine that the calculations in D = 2 + 1 would be easier 
than in D = 3 + 1 this is not really so. For example, one requirement that a first order 
transition should be well defined in a practical calculation is that the rate of tunnelling should 
be small enough for the co-exisitence of the ordered and disordered phases to be readily visible. 
This rate is governed by a factor exp{— 2<j cd L® ~ 2 /T c } = exp{— 2a°(L s / L t ) D ~ 2 } where a c d is 
the interface tension and we have written o C( i = <j q T//'~ 1 . This factor gives the suppression of 
those intermediate configurations in a tunnelling where the volume is split into two co-existing 
phases by domain walls that extend right across the lattice. So all other things being equal, 
i.e. if the dimensionless tension <t° is the same, then the size L s /L t in 2 space dimensions 
should equal the square of the size, (L s /L t ) 2 , in 3 space dimensions in order to achieve the 
same tunnelling suppression. So the typical factor L s /Lt ~ 3 in D = 3 + 1 needs to become 
a much more challenging L s /L t ~ 10 in D = 2 + 1. For reasons like this, the accuracy and 
range of our calculations is not much more impressive than in D = 3 + 1 [2] . 

4 Results 

The basic step is to calculate the maximum of x/V-> or C((3), for some given values of N, 
L t and L s . To do so we typically evaluate x/V a ^ 4 or more values of (3. These values are 
judiciously chosen so as to straddle the maximum and to be close enough to that maximum for 
each Monte Carlo run to have a usefully large number of tunnelling events between confined 
and deconfined phases. (All this clearly requires some preparatory work.) To obtain the 
maximum accurately we need a smooth interpolation between the calculated values, and this 
is provided by standard reweighting techniques (see [12] and Appendix A of [8]). In Fig. [T]we 
show an example of such a reweighted curve, for a 25 2 3 lattice in SU(6). 

We now present our results. We begin with SU(2) and SU(3) where the deconfinement 
transitiona are well known to be second order [3]. We then deal with SU(iV > 5) where they 
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turn out to be clearly first order. Finally we dwell in more detail on the delicate case of 
SU(4). In each case we perform calculations for L t = 2, 3, 4, 5 and from these we obtain the 
continuum limit of T c /y/a and (where relevant) Lh/T^. The continuum limit can be obtained 
from a standard weak coupling extrapolation of the form 

aT c (a) T c (a) T c (0) ,2 2,4 4, / 1fi x 
/ 7 \ = / / x = / TTvT + cia fji +c 2 a fi +■■■ . 18 

Here \i is some convenient physical energy scale, such as ■sja or T c , in which case a 2 T c 2 = l/I/ 2 . 
The string tension is obtained by interpolation of the results presented in [TJ. In the final 
subsection we examine the iV-dependence of these quantities as well as the scaling of finite 
volume corrections. 

One might reasonably worry that at L t = 2 the lattice spacing is too coarse for the weak 
coupling expansion in eqn ffTH]) to be applicable. (In units of T c it will be a = 0.5T C .) In fact 
there is a strong-to-weak coupling crossover in D = 2 + 1 SU(iV) gauge theories that becomes 
a phase transition at N = oo [32] in a manner very similar to the Gross- Witten transition in 
D = 1 + 1 [13]. This transition occurs at a value of (3 that is about 10% below (3 c (L t = 2). Thus 
it is not excluded that at this j3 a weak-coupling expansion is applicable. We also note that 
for light glueball masses the ratio mc/y/u appears to be described by such a weak-coupling 
expansion at comparable values of f3 [TUl [15] . For these reasons we shall (cautiously) attempt 
continuum extrapolations from L t = 2. 



4.1 SU(2) and SU(3) 

It is known from earlier calculations [3] that both SU(2) and SU(3) have second order decon- 
fining transitions, and there is good evidence that these are in the universality classes of the 
2D q = 2 and q = 3 state Potts models respectively, as one would expect from the general 
arguments of [16]. Our calculations corroborate the second order nature of the transitions. 
As an example we show in Fig. [2] a histogram of values of l p as obtained in SU(2) on a 90 2 3 
lattice at (3 = 4.9, which is very close to j3 c (L t = 3). We observe a single broad distribution 
which shows no sign of the two peak structure that would indicate co-existing phases. 

For each value of L t = l/aT c in the range L t e [2, 5] we calculate /3 C (V) for several values of 
V and extrapolate to V = oo using eqn f|T2l with v = 1,5/6 for N = 2 and N = 3 respectively. 
(These being the exponents of the relevant Potts models p2].) We show an example of such 
an extrapolation in Fig. O it is for SU(2) and L t = 3. In Table [1] we list the values of 
(3 c (y = oo) for SU(2) that have been obtained in this way, and the range of volumes used in 
each extrapolation. We also give the value of the string tension [I] that corresponds to each 
(3 C {V = oo) for each L t . (The error on a^/a includes the error on /3 C .) From this we obtain 
the values of T c j^/a = \ja^foL t shown in Table [TJ In Table [2] we display the corresponding 
results for SU(3). 
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4.2 SU(iV > 5) 

For iV = 5 it is known that the transition is first order [3, |5] and our calculations show that 
this is the case up to N = 8. We find that the strength of the transition is 'normal' and that 
the latent heat is approximately constant for N > 5. This provides convincing evidence for 
the claim that the transition remains first order for larger N and, in particular, for N = oo. 

To illustrate the first order nature of the transition we show in Fig. [4] a histogram of values 
of lp as obtained in SU(8) on a 20 2 3 lattice at value of (3 which is very close to /3 c (L t = 3). We 
observe a pronounced two peak structure that clearly indicates the co-existence of ordered and 
disordered phases. Moreover one finds that this two-peak structure becomes more pronounced 
with increasing V. 

For each value of L t = l/aT c in the range L t £ [2, 5] we calculate (3 C {V) for several values 
of V and extrapolate to V — oo using eqnflS]). We show an example of such an extrapolation 
in Fig. \5\ it is for SU(6) and L t = 3. In Tables IH 0, El we list the values of (3 C (V = oo) that 
have been obtained in this way, and the ranges of volumes used in each extrapolation. We 
also give the values of the string tension pQ that correspond to each f3 c {V = oo) for each L t . 
From this we obtain the values of T c j ^Ja = \ja^foL t shown in the Tables. 

In Table [7] we list our fitted value of the coefficient h of the 0(1 /V) correction to (3 C (V) in 
eqn(ED. 

To obtain the specific heat, Lh, as given by eqn ffTTl) . we calculate the quantity A Up defined 
in eqn( fl6l) . For small L t we can obtain it from the specific heat peak as in eqn([7]). For 
larger L t this becomes inaccurate and we obtain it by identifying sequences of confined and 
deconfined lattice fields, as described at the end of Section [3) Our results are listed in Table [HJ 
Values obtained from the specific heat are labelled by a star. In most cases where we have 
values from both methods, they are reasonably consistent. In one case this is clearly not so: 
this serves as a warning that, at least for this rather delicate quantity, our claimed control of 
systematic errors is likely to be optimistic. 

4.3 SU(4) 

If the SU(4) transition were second order, one would expect [HUH] it to be in the universality 
class of the symmetric Ashkin- Teller model (see [18] for a recent review) of which the q = 4 
Potts model is a special case. This has continuously varying critical exponents and logarithmic 
effects which, as emphasised in [4], creates difficulties in the analysis. The results of [4], 
favoured the Potts model class, although their L t = 2 transition appeared to be first order. 
Our preliminary results [5] pointed to the transition remaining first order at larger L t and this 
was confirmed by our more extensive calculations in [8]. A simultaneous recent study of this 
question presents quite convincing evidence that the transition is first order [7] . In this paper 
we present some of our results from [8] as well as some from larger volumes. 

We begin by showing the histogram of the Polyakov loop on very large volumes, at our 
smallest values of the lattice spacing. In Fig. Owe show the results obtained on a 60 2 5 lattice 
and in Fig. [7| those obtained on a 120 2 5 lattice. We see a clear 2-peak structure reflecting 
the co-existence of ordered and disordered phases. The structure becomes more pronounced 
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as L s /Lt increases from 12 to 24 and V increases by a factor of 4. The gap between the two 
maxima does not appear to shrink significantly as V increases. All this provides rather strong 
evidence that the transition is first order and remains so in the continuum limit. 

In Fig. [8] we show an extrapolation of /3 C (V) assuming an 0(1/V) correction that assumes 
it is indeed first order. As we see, such first order finite size scaling works very well. In Table [3] 
we list our extrapolated values of f3 c {V = oo), and the corresponding values of T c j^/o. In 
Table[7Jwe list the coefficient of the 0(1 /V) correction in eqnflB]). Finally in Table[8]we list the 
values of the discontinuity in the average plaquette across the transition. These will translate 
into values of the specific heat through eqn lfPTj) . 

4.4 continuum limits 

To extrapolate our calculated values of T c /^/a to the continuum limit we use eqn flTB]) with 
a/i = a^/o. We could equally well choose to use a\i = aT c = 1/L t , and we have checked that 
the results of the two types of fits are within one standard deviation of each other. (We also 
find that the quality of the fits is comparable.) In Table [9] we show the results of two such fits: 
one is an 0(a 4 ) extrapolation to L t > 2 and the second is an 0(a 2 ) extrapolation to L t > 3. 
One finds that 0(a 2 ) extrapolations to L t > 2 are statistically unacceptable. 

The best fits are rather poor, but not unacceptable, for SU(2) and SU(3) and this presum- 
ably reflects the presence of a very smooth strong-weak coupling crossover at these smaller 
values of N. For SU(5) the fit is extremely poor, and this is due entirely to our L t = 5 value. 
As we see in Table [7] the fitted coefficient of the 0(1/V) correction is anomalously large. It 
seems clear that this value happens to be out by several standard deviations. We therefore 
also show a fit (denoted by a *) that uses the values of (3 C (V = oo) given in [7J for L t = 3, 4, 5 
together with our value for L t = 2 (since this is not calculated in [7J). This gives a fit with an 
acceptable \ 2 - m f ac t our values for L t = 3,4 are consistent with those of [7J; it is only our 
L t = 5 value that is anomalous. This is reflected in the fact that the continuum extrapolation 
of our results is very similar to that obtained using the values of [7J, as we see in Table |9j 
What we can therefore say is that for larger N one has good 0(a 4 ) continuum extrapolations 
even from L t = 2. As an example, we display the SU(8) case in Fig.0 The residual systematic 
error reflected in the small difference between the 0(a 2 ) and 0(a 4 ) fits would be removed by 
calculations at slightly larger L t . 

Turning now to the latent heat, we find that we can obtain good continuum extrapolations 
from L t > 2 using just a 0(a 2 ) correction: 



We normalise by T c because that is a characteristic energy scale of the transition, and by 
iV 2 because one expects Lh oc iV 2 at large N. We choose to take the third root so that 
we extrapolate a dimensionless ratio that is of the form [m]/[m]. The results are shown in 



The values of the coefficient h at larger values of L t are too inaccurate, as we see from 
Table [7J to allow an attempt at a useful continuum extrapolation of that quantity. 




(19) 



Table M 



S 



4.5 N dependence 



In Fig. [10] we plot the continuum latent heat, in units of T c and normalised to N 2 , against 
1/N 2 , which is expected to be the leading correction to the N = oo limit. We obtain a good 
fit to N > 5 with 



It is not possible to include the SU(4) value. Indeed in our natural units, the SU(4) latent 
heat Lh is about one-half of its value at larger N, where the variation with iV is essentially 
negligible. That is to say, the deconfining transition in SU(4) is weakly first order. 

The result of attempting to fit our continuum values of T c j y/a with just a leading 0{1/ N 2 ) 
correction is shown in Fig. [TTJ (We use the values obtained by 0(a A ) extrapolations from 
Lt > 2.) It is remarkable that one obtains a good fit 



to all values of N > 2 despite the fact that the N = 2, 3 transitions are second order and the 
remainder are first order. 

It is interesting to see if the finite volume effects disappear with N 2 , as one expects if the 
free energies and interface tensions scale as oc iV 2 . A relevant quantitiy is the parameter h 
that is the coefficient of the 0(1/V) correction to the ^/-dependence of T c . As we see from 
Table [7] the errors are too large at larger L t to allow a continuum extrapolation. Instead we 
take the values at L t = 3, i.e. aT = 1/3T C , where the errors are still small and a is perhaps 
not too coarse. Plotting the corresponding values of hN 2 in Fig. [12] we see very nice evidence 
that finite volume corrections vanish as h oc 1/N 2 just as counting arguments suggest. 

5 Conclusions 

We have found that, just as in D = 3 + 1 [2], SU(iV) gauge theories in D = 2 + 1 have 
a deconfining transition that is second order at low N and first order at larger N. Again 
as in D = 3 + 1 [2], one finds that the variation with N of T c /y/a is small and that, quite 
remarkably, one can fit this ratio for all iV > 2 with just a leading 0(l/iV 2 ) correction with a 
modest coefficient, as we see in Fig. [11] 

We have seen that the SU(4) transition is almost certainly first order, although weakly so, 
as we see from its suppressed latent heat in Fig. [TO] This is reminiscent of D = 3 + 1 where 
it is the SU(3) transition that is weakly first order. 

When expressed in a natural normalisation, the value of the latent heat at larger N is 
L h ~ N 2 x (0.46T C ) 3 . By contrast, in D = 3+1 it is considerably larger [2]: L h ~ N 2 x (0.8T C ) 4 . 
This is perhaps natural given that the thermal gluons at high T have more degrees of freedom 
in three space dimensions than in two. But one should also note that if we express Lh in units 
of y/a rather than T c , most of this difference goes away. 




(20) 



T 

- 1 c 



0.9026(23) + 



0.880(43) 
iV 2 



(21) 
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The fact that the value of T c / ' yja is ~ 1.5 times larger in D = 2+1 than in D = 3+1 is quite 
natural, and already occurs for the simple Nambu-Goto string model where the deconfming 
(Hagedorn) transition occurs at T c 2 = 3/tt(D — 2) [H H9) . This simply reflects the differing 
entropy of strings or flux tubes in 2 and 3 spatial dimensions. 

It is obviously interesting to ask if the gluon plasma immediately above the phase tran- 
sition is strongly coupled, as it has been found to be in D = 3 + 1 both experimentally 
and theoretically. (See [20] and references therein.) Here the effective dimensionless 't Hooft 
coupling on the scale T c will be 

9 2 N g 2 N 

t- ~ 5 (22) 

J-c V a 

using the values in pQ. This is large enough to raise the possibility of a 'strongly coupled 
gluon plasma' in D = 2 + 1 and to motivate a numerical exploration, but that goes beyond 
the scope of the present paper. 
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SU(2) 


aT c 


L s e 


Pc 


a^Jo 


T e /y/a 


1 

2 


[15,60] 


3.4475(36) 


0.4902(13) 


1.0200(27) 


1 
3 


[22,90] 


4.943(13) 


0.3164(10) 


1.0535(33) 


1 
4 


[30,80] 


6.483(26) 


0.2317(11) 


1.0790(51) 


1 
5 


[50,85] 


8.143(57) 


0.1799(14) 


1.1117(87) 



Table 1: Critical (3 at various values of L t = l/aT c , with resulting values of T c j^/o. String 
tensions from pp. All for SU(2). 



SU(3) 


aT c 


L s e 






T /y/a 


i 

2 


[15,60] 


8.1489(31) 


0.5641(17) 


0.8864(27) 


1 

3 


[23,90] 


11.3711(91) 


0.35730(54) 


0.9329(14) 


1 

4 


[30,120] 


14.717(17) 


0.26098(41) 


0.9579(15) 


1 

5 


[40,70] 


18.131(61) 


0.20517(79) 


0.9748(38) 



Table 2: As in Tabled] but for SU(3). 



SU(4) 


aT c 


L s e 








i 

2 


[13,30] 


14.8403(26) 


0.5840(38) 


0.8562(56) 


1 
3 


[20,45] 


20.377(11) 


0.37254(52) 


0.8948(13) 


1 
4 


[30,60] 


26.228(75) 


0.27195(95) 


0.9193(32) 


1 
5 


[37,60] 


32.154(71) 


0.21423(56) 


0.9336(25) 



Table 3: As in Tabled] but for SU(4). 



12 



SU(5) 


aT c 


L s e 




aJo 

V 




1 

2 


[12,20] 


23.5315(58) 


0.5973(19) 


0.8371(27) 


1 
3 


[18,25] 


32.065(21) 


0.37912(58) 


0.8792(14) 


1 
3 




32.0765(54)* 


0.37888(52) 


0.8798(12) 


1 
4 


[24,36] 


41.057(36) 


0.27729(40) 


0.9016(13) 


1 
4 




41.113(12)* 


0.27682(30) 


0.9031(10) 


1 
5 


[30,40] 


50.67(16) 


0.21626(80) 


0.9248(34) 


1 
5 




50.275(20)* 


0.21823(22) 


0.9165(10) 



Table 4: As in Tabled] but for SU(5). 



SU(6) 


aT c 


L s e 








i 

2 


[10,20] 


34.1848(52) 


0.6033(35) 


0.8288(48) 


1 
3 


[15,25] 


46.395(17) 


0.38058(94) 


0.8759(22) 


1 
4 


[20,30] 


59.372(48) 


0.27883(46) 


0.8966(15) 


1 
5 


[30,37] 


72.59(24) 


0.21975(88) 


0.9101(36) 



Table 5: As in Tabled] but for SU(6). 



SU(8) 


aT c 


L s e 




a^Jo 


T c /y/a 


i 

2 


[5,14] 


61.245(11) 


0.60688(92) 


0.8239(13) 


1 
3 


[10,20] 


82.704(43) 


0.38525(54) 


0.8652(12) 


1 
4 


[10,18] 


105.52(11) 


0.28190(43) 


0.8868(14) 


1 
5 


[12,25] 


128.30(28) 


0.22332(61) 


0.8956(25) 



Table 6: As in Tabled] but for SU(8). 
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h 


aT c 


SU(4) 


SU(5) 


SU(6) 


SU(8) 


1 

2 


0.275(14) 


0.190(15) 


0.132(7) 


0.0567(20) 


1 

3 


0.390(39) 


0.239(36) 


0.188(13) 


0.079(11) 


1 

4 


0.50(19) 


0.212(50) 


0.204(29) 


0.066(16) 


1 

5 


0.88(24) 


0.71(15) 


0.25(14) 


0.052(29) 



Table 7: The coefficient of the 0(1/V) term in the variation of /3 C (V) in eqn(|HJ). 



A Up at T = T c 




SU(4) 


SU(5) 


SU(6) 


SU(8) 


2 


0.00603(9)* 
0.00631(37) t 


0.00886(11)* 


0.01118(10)* 


0.01238(9)* 


3 


0.00073(8)*+ 
0.00087(6) 


0.00117(12)* 
0.00125(4)+ 


0.00138(9)*+ 


0.00206(6)* 
0.00146(3)+ 


4 


0.00024(5) 


0.000336(9) 


0.000352(20) 


0.000377(24) 


5 


0.00010(1) 


0.000140(5) 


0.000148(4) 


0.000157(8) 



Table 8: The difference between the average plaquette in the confined and deconfined phases 
at T = T c . Entries with a * are calculated from the specific heat. Where there are two values, 
the f indicates the value used in the continuum fit of the latent heat. 



continuum limit 






0(a 4 ) fit L 


t >2 




0{a 2 ) 


fit L t > 3 




N 




Cl 


c 2 


X 2 /n d f 




Cl 


X 2 /n d f 


2 


1.1238(88) 


-0.87(14) 


1.83(45) 


3.3 


1.1224(90) 


-0.70(15) 


3.3 


3 


0.9994(40) 


-0.65(6) 


0.94(15) 


2.1 


0.9890(31) 


-0.44(4) 


1.9 


4 


0.9572(39) 


-0.56(5) 


0.76(13) 


0.1 


0.9516(33) 


-0.41(3) 


0.5 


5 


0.9404(38) 


-0.53(5) 


0.68(11) 


8.2 


0.9325(28) 


-0.38(3) 


11.1 


5* 


0.9380(19) 


-0.49(3) 


0.58(8) 


0.9 


0.9337(14) 


-0.38(2) 


3.4 


6 


0.9300(48) 


-0.46(7) 


0.51(16) 


1.1 


0.9228(35) 


-0.33(4) 


1.0 


8 


0.9144(41) 


-0.38(5) 


0.38(11) 


0.4 


0.9114(25) 


-0.31(3) 


0.1 



Table 9: Continuum limits of T c /yja for various N, using either a 0(a 4 ) or 0(a 2 ) formula, 
with values of the coefficients of the lattice corrections, and the quality of fit. 
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latent heat 


N 


i (L h \i 

T c \WJ 


X 2 /n df 


4 


0.393(13) 


1.4 


5 


0.448(4) 


1.4 


6 


0.459(4) 


0.9 


8 


0.453(4) 


0.9 



Table 10: Continuum limit of the latent heat using an 0(a 2 ) extrapolation to 2 < L t < 5. 



Tj^a vs TV 


0(l/N 4 ) fit 


0(l/iV 2 ) fit 


N > 


N = oo ci c 2 X 2 / n d/ 


N = oo ci X 2 /n df 


a 4 2 


0.9035(41) 0.85(11) 0.12(31) 0.23 


0.9026(23) 0.88(5) 0.18 


a 2 2 
a 2 3 


0.9016(26) 0.76(7) 0.46(33) 0.70 
0.8977(38) 0.94(14) -1.0(1.1) 0.07 


0.8989(17) 0.85(4) 0.98 
0.9007(20) 0.81(5) 0.35 



Table 11: Results of fitting the A^-dependence of T c /^a either to 0(1/ TV 4 ) or to 0(1/A^ 2 ). 
Results shown separately for the different extrapolations to the continuum limit. 
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Figure 2: Number of lattice fields, N, with a given value of the Polyakov loop \l p \. In SU(2), 
at = 4.9, on a 90 2 3 lattice. 



4.95 




0.05 0.1 0.15 

(L s /L t )-^ 



0.2 



Figure 3: A plot of the critical coupling f3 c (V) against (L s /L t ) x l v with an extrapolation to 
V = oo. In SU(2) for aT c = 1/L t = 1/3. 
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0.5 



1.5 

\l P \ 



2.5 



Figure 4: Number of lattice fields, N, with a given value of the Polyakov loop \l p \. In SU(8), 
at = 82.5, on a 20 2 3 lattice. 



46.4 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 

{L t /L s f 



Figure 5: A plot of the critical coupling y = (3 C (V) against x = (L s /L t ) 2 with an extrapolation 
to V = oo. In SU(6) for aT c = 1/L t = 1/3. 
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i = 31.78 

13 = 31.9 

(3 = 32.0 




0.35 



WIN 



Figure 6: Histogram of the number of lattice fields with a given value of the Polyakov loop 
\l p \/N. In SU(4), at the (3 indicated, on a 60 2 5 lattice. 




Figure 7: Histogram of the number of lattice fields with a given value of the Polyakov loop 
|Z p |/iV. In SU(4), at the (3 indicated, on a 120 2 5 lattice. 
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0.05 



Figure 8: A plot of the critical coupling y = f3 c (V) against x = (L s /L t ) 2 with an extrapolation 
to V = oo. In SU(4) for aT c = 1/L t = 1/3. 
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0.8 



0.6 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

0?(J 



Figure 9: Values in SU(8) of T c (a)/y/a(a) plotted against a 2 a and extrapolated to the con- 
tinuum limit using either linear (...) or quadratic (-) extrapolations. 
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J I 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 

1/N 2 

Figure 10: Latent heat, with a natural normalisation, obtained at a fixed lattice spacing 
a = 1/3T C . The fit to large N uses a leading 0{1/N 2 ) correction. 

1.5 I 1 1 1 1 1 

Is. 




0.5 1 1 1 1 1 1 1 

0.05 0.1 0.15 0.2 0.25 0.3 

1/N 2 

Figure 11: Continuum values of T c /^/a plotted against 1/N 2 . The fit to all N uses just a 
leading 0(1 /N 2 ) correction. 
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Figure 12: iV-dependence of the coefficient h of the 0(1 /V) correction to (3 C (V), obtained at 
a fixed lattice spacing a = 1/3T C . 
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